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Abstract. We provide an analytical proof of the existence of a stable periodic 
orbit contained in the region of coexistence of the three species of a tritrophic 
chain. The method used consists in analyzing a triple Hopf bifurcation. For 
some values of the parameters three limit cycles bear via this bifurcation. One 
is contained in the plane where the top-predator is absent. Another one is not 
contained in the domain of interest where all variables are positive. The third 
one is contained where the three species coexist. The techniques for proving 
these results have been introduced in previous articles (see [2] 1181 ) and are 
based on the averaging theory of second— order. Existence of this triple Hopf 
bifurcation has been previously discovered numerically in |14) . 

1. Introduction 

During these last 80 years, after the seminal works of Lotka [17] and Volterra [27] . 
one of the main topics in mathematical ecology has been the study of (di) trophic 
food chains. This has been made by analyzing many different planar differential 
systems under the common name of prey-predator models, for instance see [TJ 
[7]. The existence of limit cycles, attractors, and several kind of bifurcations are 
the characteristics of those models which have been used to explain the complex 
behaviors observed in such systems. 

In the late seventies some interest in the mathematics of tritrophic food chain 
models (composed of prey, predator, and top-predator) appeared, see for example 
[U [lOl [8] and Predator-Prey models with parasitic infection QT] . The model we 
analyze in this article describes a tritrophic food chain composed of a logistic prey 
(x), a Holling type II predator (y), and a Holling type II top-predator (z). After a 
rescaling of the variables, it is given by the following system of ordinary differential 
equations (see OH E2 Q3 Q3[ [20] for more details): 

x _ aiy 
k 6i+a 



\bi + x b 2 + y I 

0.2V 



b 2 + y 



-do). 
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In order to preserve the biological meaning of the model, the 8 parameters of this 
system are assumed to be strictly positive. Similar types of systems have been 
studied in the case when time scales of the variables are different so that methods 
of approximations of slow- fast systems can be applied [25l [26]. We emphasize 
that we do need here this type of approximation. To clarify possible applications 
of our result, we focus in a range where the population of the superpredator z 
remains small (near z = 0). So in our view here the superpredator remains of small 
amplitude compared with the two other populations. Within this range, we look 
for the possibility of a periodic rhythms which is stable and where the three species 
coexist. 

The bifurcation analysis can be carried out with respect to the two parameters 
(di, d 2 ). The normal form analysis around one of the stationary points (called latter 
p 3 in this article was first carried out by Klebanoff and Hastings ([13]) and latter 
improved by [15] . These authors (Kuznetsov and Rinaldi) also discovered numeri- 
cally the existence of a strange attractor. In the limit of slow- fast systems Muratori 
and Rinaldi ([22]) developed a singular perturbation approach which supports the 
existence of a homoclinic intersection of the stable and unstable manifolds of the 
limit cycle which is contained in the (x, ?/)-plane. Classical arguments in the theory 
of dynamical systems can then be used to deduce the existence of infinitely many 
limit cycles of increasing period (but nothing is known in general on their stability) . 
See also f[25] and [26]). 

We are interested in the limit cycles of system (p} , mainly in the ones which come 
from a Hopf bifurcation. There are several papers dedicated to these limit cycles 
see for instance [HI [3j [H El El EU [16] . But in all these papers the existence of a 
triple Hopf bifurcation was not proved analytically, the results there are essentially 
numerical. The existence of the triple Hopf bifurcation was discovered in [14] . We 
prove here that there are systems ([lj having 1 or 3 small amplitude limit cycles 
coming from a Hopf bifurcation, and we show how to study the type of stability 
of such limit cycles. The tool for obtaining these results is the theory of averaging 
of second order. In fact as we shall show from the singular point bifurcate 3 small 
amplitude limit cycles one of them is contained in z > 0, the other in z = and 
the last in z < 0. Of course the unique small amplitude limit cycles which have 
biological meaning are the ones contained in z > 0. 

As a guide to the potential "ecological user" , we should emphasize that we can 
prove the existence of a Hopf bifurcation of second order provided that the logistic 
growth of the prey satisfies the two relations: 

_ bi(ai + di) 
^ (ai — d\)k ' 

and 

(ai — d 1 ) 2 (a 1 b 1 — 2b 2 d{) ' 
and that the rate of exponential decay of the superpredator satifics: 

, _ a\a 2 b\ 

2 ~ a\b 2 k + b 2 d\k + ai{b\ - 2b 2 dxk) ' 

There is also a fourth condition which defines a full open set in the parameter space: 

ai&i > 2b 2 d\. 
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These four conditions ensure the existence of a periodic orbit contained in the 
domain of coexistence of the three species. A last inequality is required to ensure 
the stability of this periodic orbit. Altough from the viewpoint of applications 
imposing inequalities on the parameters is not constraining too much (provided 
they are compatible and we check that the set of existence of a stable periodic orbit 
is not empty by giving one example), it is of course much less natural to impose 
the three first conditions above. These conditions seem necessary for the analytical 
proof. Practically it might be enough to check that the values of the parameters 
remain close to those prescribed above but this last point is beyond the dynamical 
systems techniques we use here. 

The paper is structured as follows. In the appendix [4] we summarize the basic 
results on the averaging theory of second order that we shall need for studying the 
Hopf bifurcation. The explicit results on the Hopf bifurcation are stated in Theorem 
[T] at the end of section [21 In this section is also proved Theorem [TJ An example 
showing the existence of three stable small amplitude limit cycles bifurcating from 
a singular point of system ([1]) is given in Section 3. 

2. The Hopf bifurcation 

We separate the study of the Hopf bifurcation in different subsections. 

2.1. The singular point which will exhibit the Hopf bifurcation. The dif- 
ferential system ([1]) can have the following six singular points: 

Pi = (0,0,0), 



where 



P2 = (fcp,0,0), 

&idi 6i(&idi + (di — ai)k p) 



P3 



PA = 0, 



di — di (ai — di) 2 fc 

b 2 d 2 b 2 d x 



P5 



P6 = 



a 2 — d 2 ' a 2 ~ d 2 
A + VB b 2 d 2 b 2 (ai-di)VB-b 2 C(a 2 -d 2 



2(a 2 - d 2 ) ' a 2 - d 2 ' (a 2 -d 2 )(VB + £>) 
A-VB b 2 d 2 6 2 (ai-di)V5 + 6 2 C*(a 2 -d 2 ) 



2(a 2 -d 2 )'a 2 -d 2 ' (a 2 ~d 2 )\WB-D 



A = —a 2 b\ +b\d 2 + a 2 k p — d 2 k p, 

B = 4(a 2 - d 2 )(a 2 (6i + k pf - d 2 {b\ + \a x b 2 k + 2b x k p + k 2 p 2 )), 
C = a\bi + b\d\ — a\k p + d\k p, 
D = a 2 b\ —b\d 2 + a 2 k p — d 2 k p. 

Of course such singular points exist always that their denominators are nonzero 
and the expression B > 0. 

The study of the stability of p§ and pe looks very tedious in full generality due 
to the very long expressions of their eigenvalues. So the analysis of their possible 
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Hopf bifurcations is out of scope. Since the eigenvalues of p\ and P2 are always 
real, we must put our interest in the points P3 and p^. We choose here to study P3 
because it seems more meaningful to proceed in the neighborhood of z — if we 
keep in mind the situation of biological interest we presented in the introduction. 
It could be quite possible to proceed with p^. 

We shall study the Hopf bifurcation at the singular point p%. The eigenvalues at 
this singular point are 

\±= - — -, — rr f — a\b\d\ — bid 2 + a\d\k p — d\k p ± \fK 

2ai(ai — di)k V 

a 2 &i(6idi + (— Oi + di)k p) 

P — — «2 1 



b\di — 62(01 — di) 2 k + b\ (— ai + d\)k p 
where 

A = di(-4ai(ai - e?i) 2 fc(-Mi + (a x - d x )k p) + di(ai(&i - k p) + d 1 (b 1 + kp)f). 

It is well known that a necessary condition in order to have a Hopf bifurcation 
at j»3 is that their pair of complex eigenvalues A± when A < must cross the 
imaginary axis. So we will take 

-aibidi — bidf + aidik p — dlk p 2 

(2) Ke(A±j = — TT7 =e I, 

la\{ai — d\)k 

where e is a small parameter necessary for applying the averaging theory in order to 
study the Hopf bifurcation and I is an arbitrary parameter. Initially we had taken 
instead of e 2 l the expression el\ + e 2 h- But applying the averaging theory of first 
order we do not obtain any information about the Hopf bifurcation and in order 
that the averaged function of first order _Fio (see the appendix) becomes identically 
zero we must take l\ = 0. This is the reason that now we are taking directly e 2 l 
instead of el\ + e 2 l 2 - 

In order to apply the averaging theory we must write our differential system ([1]) 
into the normal form of the averaging, i.e. into the form (| 1 2[) . Then we need to 
take 

(3) /i = e 2 m, 

where m is an arbitrary parameter. As in the equation ((2]) at the beginning we 
took £mi + £ 2 iri2 , but for the same reason than before we must take mi = 0. In 
short we are taking directly e 2 m instead of emi + e 2 m,2- 

Solving equalities and ([3]) with respect to c? 2 and p we obtain that 

a\a2b\d\ + Ee 2 + 2a\bi{di — ai)klme A 



d2 = 



(4) 



1 d\ (b 2 kaj + (b 2 - 2b 2 d 1 k) a x + b 2 d\k + 2ai&i(ai - dx)kle 2 ) ' 
2a\{a\ — d\)kle 2 + bidi(ai + d\) 



(ai — di)d\k 

with E = —b2kmd\ + a\ \—mb\ — 2a2klb\ + 2b2d\kvn)j d\ + a 2 k(2a,2bil — b2dim). 
Now the eigenvalues of the singular point p^ become 



2] Jk (e 2 k I (le 2 - 2a t + 2d!) - Ml) ? 
e I ± — ; , e m. 
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2.2. The system in the normal form for applying the averaging theory. 

One of the most difficult steps for applying the averaging theory to a given differ- 
ential system, as our system ([1} in order to study its limit cycles, is to write it in 
the normal form (|12p of the averaging theory. Moreover in our case we want to 
study the small amplitude limit cycles coming from a Hopf bifurcation. For doing 
all this we shall follow the next steps: 

(i) First we translate the singular point P3 at the origin of coordinates doing the 
change of variables (X, Y, Z) = (x,y, z) — p^. 

(ii) We write the linear part at the origin of the differential system (X,Y,Z) in its 
real Jordan normal form. For this we do the change of variables 




/ dik 

1 — r 1 

e 2 kl ^ e 2 (2l-m) 



\ 



di 
G 




where 



F : 
G- 
H 



y/hdi. - k (e 2 k I (le 2 - 2a x + 2di)), 

HI/ [aia 2 bidik (2(oi - d x )kle 2 + Mi)] , 

b 2 kd{ + oi(6f - 2e 2 klb 1 - 2b 2 d 1 k)d 1 + a\k{2b x le 2 + Mi) 



I = k (m(m - 2l)e 2 + 2a x l - 2d x l) 



Ml 



(iii) We write the new differential system (U, V, W) in the cylindrical coordinates 
(R, 9, W) defined through U = Rcos 9, V = i?sin 6 and W = W. 

(iv) In order to study the small amplitude limit cycles around the origin of coordi- 
nates we do the rescaling (R, 9, W) = (er, 6, ew). 

(v) Finally the system (f, 9, w) is written as (r', w') = [dr/dO, dw/d6) and we obtain 
the tritrophic food chain model in the normal form of averaging: 

r' = sF n (9,r,w) + e 2 F 21 (9,r,w) + 0(e 3 ), 
w' = eF 12 (9,r,w) + e 2 F 22 (9,r,w) + Q(e 3 ), 



where 



-Fii — Ri/Tq, 

F 2 \ = (R 2 Tq — R±Ti)/T , 

F 21 = Wx/Tq, 

F 22 = (W 2 T - WxT^/T 2 , 
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and 



i?i = ^ cos 3 - , cos 2 sin =— cos 2 0+ 

awfbidik bi a\k 

^h^L cos e + ,(^5+ 

aifc 2 cuVhdJk (b 2 ka 2 + b\ ai - 2b 2 d 1 ka 1 + b 2 djkf 

a 2 b 1 b 2 k 2 a\ - Zb\d x k 2 a\ + 26 2 fe 2 fc a 4 + fcfaf + 2b 2 d\k 2 a\ - 3a 2 &iMifc 2 a?- 
2b\b 2 d 1 ka\ + 2b\d\k 2 a\ + Za 2 b 1 b 2 d\k 2 a\ + fcfdia 2 - 2b\b 2 d\ka\- 
ibld\k 2 a x - a 2 6i6 2 difc 2 ai + 2b\b 2 d\ka 1 + b\d\ k 2 )rw cos sin 0+ 
(ai — — a^fcr 2 



a\bi^/b\dik 



■ cos sin 0+ 



■sin sin ( 



6 2 fca 2 + 6 2 ai — 2b 2 d\ka\ + b 2 d\k a\k 



b 2 k 2 r 2 w{ai - dif . 3 d ir 3 ( ai - dx) 2 4 

R 2 = ~ sm J — cos 4 0+ 

a 1 b 1 {b 2 ka 2 + b 2 a 1 -2b 2 d 1 ka 1 + b 2 d\k) a{b{ 

w 3 (ai-di) 2 . r(w 2 a5 - 3di«; 2 af - 6 2 Za 2 + 3d 2 w 2 ai - d?w 2 ) . 2 

(ai - di)di(di - ai)kr 3 , . 3(ai - di) 2 \/bidikwr 2 , 

„ , — cos 3 sm + — Vpr- COS 0- 

aibf^bidik afbfk 

(ai- dxfdxkr 3 2 . 2 (on - di) 2 (2ai + (ii)rar 2 2 . . 

^ cos 2 sin + „,„ cos 2 sm 0+ 

a 2 &j a 2 6 2 

(— 3w 2 a? + bikla 2 + 6diw 2 ai — 3o!?w 2 )r „ J 2 
T — cos z + — cos sm 0+ 



T Q = - 



^bidikw 3 (ai - di) 2 (ai + 2cii)ru; 2 (ai - di) 2 

97 , , „ — cos 5 . cos sm 0, 

a\bxd\k l a\b\\Jb\d\k 

\Jb\d\k 
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(ai - di)dir o 2(dx - a^^Jbidikw 2 

Ti = : cos J 9 — COS ( 

a\bi a\bik 



Jb\dik(ka 2 — 2d\ka\ + b\d\ + d 2 k)r „ . . {a\ — d\)w 

-j- cos z 9 sin 9 ; cos 9+ 

aib^k a\kr 

{b\ka\ — 2b\d\ka\ + 2b 2 d\ka 2 + 26 2 e?iai + bid 2 kai — 4b 2 d\ka\ + 2& 2 df fc)u; 
ai£;(& 2 fca 2 + b\a\ — 26 2 <ii/cai + 6 2 G? 2 fc) 

(ai 'f l)2r cos 9 sin 2 9 - ± sin 2 9 sin 9, 

a2hb 2 (a 1 - di) 3 (ai + dx)kw . 

Wi= — 2 Sm ^' 

(& 2 fca 2 + b\a\ — 2b 2 d 1 kai + 6 2 (i 2 fc) 



cos 9 sin 9+ 



a 2 blb 2 (ai - di) 4 (a x + d x ) 2 kw . 2 

W 2 = raw - ,, , 9 p5 — — — : — pr-^ sin 9. 

{b 2 ka\ + b\a x - 2b 2 d 1 ka 1 + b 2 d\k) 3 

where 

J = k(2b 3 k 3 d\ - \2aib 3 2 k 3 d\ + i0a 2 b 3 2 k 3 d\ + 6ai6f&|fc 2 df - 40a?6|fc 3 df- 

2Aa 2 b\b 2 2 k 2 d\ + ?>Qa\b 3 k 3 d 3 + 36a?& 2 6f.fc 2 d? + a 2 a 2 b\b 2 k 2 d\ + %a\b\ b 2 kd\- 
\2a\ b 3 2 k 3 d\ - 2\a\b\b\k 2 d\ - 3a?a 2 6?6 2 A: 2 d 2 - I2a\b\b 2 kdl + 2a?6f di + 
2a\b 3 2 k 3 d 1 + §a\b 2 b 2 2 k 2 d x + 3af a 2 b\b 2 k 2 d 1 + 6afbfb 2 kd! - a\a 2 b\b 2 k 2 )r 2 w 
(ai - di) 3 , 



K = albl^bidik(b 2 kal + b\a\ - 2b 2 dikai + b 2 d\k) 3 , 

L = (ai - rfi) 2 \/6irfifc(&ifc 2 af + a 2 bib 2 k 2 a\ - ^b\d Y k 2 a\ + 2b\b 2 ka\ + bfaf+ 
§b 2 2 d\k 2 a\ - 2a 2 bib 2 dik 2 a\ - ^b 2 b 2 d x ka\ - \b 2 2 d\k 2 a\ + a 2 b 1 b 2 d 2 k 2 
2b 2 b 2 d 2 k ai + b\d\ fc 2 )w, 

M = aibidik(b 2 ka 2 + b\a x — 2b 2 d\kai + b 2 dfk) 2 . 

2.3. The computation of the small amplitude limit cycles. Now we compute 
the function Fiq, see Theorem and we obtain 

r 2lT -1 1-2-K 



F y Jr.,r)= [i-jf F u (O,r,w)d0,^f F 21 (9,r,w)d9 



Nw\ 

0.--P-). 

where 

AT = 6ifca? - 2b\dika\ - 2b 2 d\ka\ - 2b\d\a\ + bidjkai + Ab 2 d\kai - 2b 2 d\k, 

P = a\b 2 ka\ + b\a\ — 2b2d\ka\ + b 2 d\k) . 
Taking 

(5) k = 2ai6?dl 



(ai - d 1 ) 2 (a 1 bi - 2b 2 d 1 ) : 
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we obtain that the averaged function of first order is identically zero. We must 
compute the averaged function of second order, for more details see the appendix. 

We note that since in the expressions of with i,j G {1,2} appears yk we 
need that 



(6) 



aih - 2b 2 di > 0. 



From (fT5|) the averaged function of second order F20 (r, w) has the two compo- 
nents: 



F 201 (r,w) 



Nr 



2(oi — di) 2 {a\bi — 2b 2 di){bia\ — bid\a\ — ^a 2 b 2 d\)w 2 - 



2y/2a\b\d x 
albld 1 (4a 1 lb' 2 + d 1 (d 1 - ai)r 2 )), 



F 202 (r, w) 



where 



a\b\ 



(ai - d 1 ) 2 {a 1 b 1 - 2b 2 d 1 )N 3 w 



(a\m 2 b\ — d>a\a 2 b 2 d\r 2 bi — 2a 2 b 2 (a\ — d\) 2 d\{a\bi — 2b 2 d\)i 



N = 



(ai - di) 2 {aibi - 2b 2 d\) 



In order to look for the small amplitude limit cycles bifurcating from the origin 
of system (fTJ) , after all the changes of coordinates that we did and according with 
Theorem [21 we must find the zeros (ro,u>o) with ro > of the system 

(7) F 20 i{r,w)=0, F 202 {r,w) = 0. 

such that the Jacobian 

/ dF 201 

(8) 



det 



dr 
dF 2 Q 2 



dF 2Q1 \ 
dw 



\ dr 



9F 2 q 2 
dw I 



(r.w) = (r ,Wo) 



be nonzero. 



It is easy to check that system ([7]) can have at most 3 solutions satisfying ([5]) 
according with the values of the parameters of system |T|). More precisely one 
solution (ri,wi) is 



(9) 



n =2b u 



1 ail 
(ai - 



2b iy / R 1 , wi=0; 
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of course it exists if R\ > 0. The other two possible solutions are (r2,±u;2) where 



ai&i I a\b\ (a\bi(ai — di)m — Aa 2 b 2 d\(l + m)) a\b\ p^- 
d\ y a 2 b 2 di {bb\a\ — hbidi<x\ — 2Aa 2 b 2 d\) d\ 2 



(10) 



U'2 



ajbf / 24a 2 b 2 d(l - albi(a! - di)m 



V2 V 0262(01 - d 1 ) 2 d 1 (2b 2 di - ai&i) (-5&iaf + db^aj + 2Aa 2 b 2 d\) 

w 2 - 



V2 



of course again these last two small amplitude limit cycles will exist if R 2 > and 
W 2 > 0. 

When the parameters of system ([1} are such that we have 3 small amplitude 
limit cycles, then the one which has initial conditions on the plane w = (see 
Theorem [5]) remains in this plane because it is invariant by the flow of the system. 
Again since this plane is invariant, the other two small amplitude limit cycles which 
have initial conditions on half-spaces w > and w > remain in such half-spaces. 
Of course the invariant plane w = corresponds to the plane z = in the initial 
coordinates. 

According with the statement (b) of Theorem [2] we can compute the type of 
stability of these small amplitude limit cycles computing the eigenvalues of the 
matrix |8]). More precisely, let Ai and A2 be the two eigenvalues of the matrix ([8|) 
evaluated on a zero (ro,u>o) of system ([7]), then the small amplitude limit cycle 
associated to the zero (r , w ) 

(I) is a local repeller if Re(Ai), Re(A2) > 0, 
(II) is a local attractor if Re(Ai), Re(A2) < 0, 

(III) has two invariant manifolds, one stable and the other unstable, which locally 
are formed by two 2-dimensional cylinders. 

We note that conditions (|3|) with e — 0, © and ^ are necessary in order to 
apply the second order averaging theory. More precisely we need that system (fTJ) 
satisfies the conditions 

a\a2b\d\ 



d 2 = 



(11) 



k = 



a\b 2 d\k + b 2 d\k + aid\(b\ — 2b 2 dik) ' 

&i(oi + di) 
(a\ — d\)k ' 

2a\b\d\ 



(ai - di) 2 (ai6i — 26 2 di) ' 
0< a 1 b 1 -2b 2 d 1 , 

in order that we can apply the averaging theory of second order for studying its 
small amplitude limit cycles. 

In short we have proved the next result. 



Theorem 1. The following statements hold. 
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(a) The tritrophic food chain model given by system ([I]) satisfying the four 
conditions (llip has a triple Hopf bifurcation at the singular point p% if R\ , 
R 2 and W 2 are positive (these last three expressions are defined in ^ and 

(Hhd;. 

(b) More precisely, under the assumptions of statement (a) three small am- 
plitude limit cycles bifurcate from P3 with initial conditions in the coor- 
dinates (U, V, W) given by {eri,Q,ew\) and {er 2 , 0, ±£^2), where the val- 
ues of Ti and Wi are given in (JSJ) and (|10p . Moreover in the variables 
{x,y,z) the small amplitude limit cycle coming from the initial conditions 
(eri, 0, ew\) lies on the plane z = 0, the one coming from the initial condi- 
tions {er 2 , 0, sw 2 ) lies in the half-space z > 0, and the remainder one lies 
in the half-space z < 0. 

(c) Under the assumptions of statement (a) the kind of stability of the small 
amplitude limit cycles is determined by the eigenvalues of the matrix ([5} as 
it is explained in (I), (II) and (III). 

3. An example 
When e = the relations (fTTj) for the values 

ai = 5, a 2 = 0.1, h = 3, 6 2 = 2, di = 0.4, 

become 

d 2 =0.09, p = 27.74, k = 0.13, aib x - 2Mi = 13.4. 

These values are compatible with the biological conditions. Moreover the parame- 
ters I and m are free. Taking / = 400 and m = lwe obtain that the eigenvalues of 
the matrix in the expression ((8|) are 

-154.96 and -0.32 for (r u wi) = (221.16, 0), 
-135.14 and -0.29 for (r 2 , w 2 ) = (207.24, 39), 
-135.14 and -0.29 for (r 2 , w 2 ) = (207.24, -39). 

Therefore the three small amplitude limit cycles are local attractors for these values 
of the parameters. 

4. The appendix: Averaging theory of second order 

We shall use the following result. 

Theorem 2 (Second order averaging method). We consider the following differ- 
ential system 

(12) ^ = eF x {t,x) +e 2 F 2 {t,x) +e 3 R(t,x,e), 

dt 

where F 1: F 2 : R x D — > R", JJ:lxDx (—£f,Ef) — > R" are continuous functions, 
T -periodic in the first variable, and D is an open subset o/R™. We assume that 
(i) Fi(t,-) e C^D) for all t e R, Fx, F 2 , R and D x Fi are locally Lipschitz 
with respect to x, and R is differentiable with respect to e. 
We define F w , F 20 : D — > R" as 

1 f T 

F w{z) = y J F 1 (s,z)ds, 
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■t r 



D z Fi{s,z)- I F 1 (t,z)dt + F 2 (s,z) ds, 



o 



and assume moreover that 
(ii) for V C D an open and bounded set and for each e € (— £/, £/) \ {0}, there 
exists a e <E V such that Fio(a £ )+eF2o(a £ ) = and ^(-Fio+e^o, V, a e ) ^ 0. 
Then the following statements hold. 

(a) For \e\ > sufficiently small there exists a T '-periodic solution Lp(t,e) of 
system (|12|) such that (/?(0, e) = a e . 

(b) // the function Fiq + ei^o * s C 1 the stability or instability of the limit cycle 
(p(t, e) is given by the stability or instability of the singular point a e of the 
averaged system 



corresponding to system (| 1 2() . In fact, the singular point a e of the averaged 
system has the stability behavior of the Poincare map associated to the limit 
cycle <j)(t, e). 

The proof of statement (a) can be found in [2J, and of statement (b) is standard, 
sec for instance I35 ) US] . 

The expression dsiFiQ + eF2o, V, a e ) ^ means that the Brouwer degree of the 
function Fiq + eFw ■ V — ► R™ at the fixed point a e is not zero. For instance it is 
sufficient that the Jacobian of the function F±o + e-Pjo a t a e is not zero in order 
that the mentioned Brouwer degree be nonzero. 

If F\o is not identically zero, then the zeros of Fio + ei^o are mainly the zeros of 
Fiq for e sufficiently small. In this case the previous result provides the averaging 
theory of first order. 

If Fio is identically zero and F20 is not identically zero, then the zeros of -Fio+ei^o 
are mainly the zeros of -F20 for e sufficiently small. In this case the previous result 
provides the averaging theory of second order. 
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